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We report a direct lattice calculation of the K to 7T7Z; decay matrix elements for both AI = 1/2 
and 3/2 channels on 2+1 flavor, domain wall fermion, 16 3 x 32 lattices with zero %% relative 
momentum and m K = 420 MeV. All K° to KK contractions are carefully listed and calculated. 
The decay into the isospin zero nn final state, which receives contributions from the disconnected 
graphs, is very difficult to calculate, but a clear signal in the similar disconnected nn correlator can 
be seen. Preliminary results, some with large errors, will be presented for the various contributions 
to the renormalized weak matrix elements Aq and Ai. We obtain Re(Ao) with 25% error in the 
case of zero momentum on shell decay, and find a factor of 6 enhancement for the AI = 1/2 rule 
in the 420 MeV pion system. 
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1. Introduction 

The observation of the AI = 1/2 rule and direct CP violation in kaon decays has attracted 
extensive theoretical study but both still lack a quantitative explanation. The running of the Wilson 
coefficients down to ~ 2 GeV which represents short distance physics can only explain a factor of 2, 
far less than the experimental factor of 25 AI = 1/2 rule enhancement. The remaining enhancement 
comes from hadronic matrix elements which requires non-perturbative treatment. On the other 
hand, direct CP violation in kaon decays serves a very important check of the standard model's 
CKM mechanism of CP violation. While experimentalists have measured Re(£ ; /£) = 1.65(26) x 
10" 3 , with only 16% error, there is no reliable theoretical calculation based on the Standard Model. 
Previous chiral perturbation theory based lattice QCD calculation using 2+1 dynamical domain 
wall fermion fails to give a conclusive result because of large systematic error [1]. Therefore, a 
direct lattice calculation of K — > nn decay is extremely important to provide an explanation of the 
nature of AI = 1/2 rule and direct CP violation. This is a notoriously difficult calculation, but with 
the increasing advance of computing power, we want to show that it is now accessible. 

In this paper, we try to do a direct, brute force calculation of the needed weak matrix elements. 
The isospin n — 71 final state involves disconnected graphs that make the calculation very difficult. 
For these graphs, the noise does not decrease with the increasing time separation of the source and 
sink, while the signal does. Therefore, huge statistics is needed to get a clear signal. As a trial 
calculation, we do this on a relatively small lattice , so it is easier to collect large statistics. We 
concentrate on the study of the statistical uncertainty since it is the major difficulty of the problem. 
We will mainly report our results for the AI = 1/2 calculation in this paper. The inclusion of the 
AI = 3/2 part is for completeness; a much better calculation of the AI = 3/2 amplitude alone on 
a large lattice can be found in [2]. In Section 2, we summarize our setup of the calculation. Then 
our 71 — 71 scattering results are given in Section 3. Section 4 shows the details of the K° to %% 
contractions, and the calculated results and conclusion are shown in Section 5. 

2. Computational Details 

Our calculation uses the Iwasaki gauge action(/3 = 2.13) and a 2+1 flavor(wi/ = 0.01, m s = 
0.032) domain wall fermion action, with space time volume 16 3 x 32, and L s = 16. The inverse 
lattice spacing for these lattices is determined to be 1.73(3)GeV, and the residual mass is m res = 
0.00308(4). The propagators are calculated on each of the 32 time slices using a Coulomb gauge 
fixed wall source (used for mesons), and a random wall source (used to calculate loops in type3 
and typeA graphs in Fig. 2). One propagator needs 12 (3 colors and 4 spins) Dirac operator 
inversions, so all together we carry out a few sets of 384 inversions for different sources and quark 
masses on a given configuration. This calculation is accelerated by a factor of 2-3 for mi = 0.01 
by computing the Dirac eigenvectors with the smallest 35 eigenvalues and limiting the conjugate 
gradient inversion to the remaining orthogonal subspace. In order to obtain an on shell K — > Tin 
decay amplitude, the valence strange quark is partially quenched and its mass is chosen to be 
m s = 0.066,0.099,0.165, with the corresponding kaon mass shown in Tab. 1. In the following 
section, we will see that we can interpolate to on shell decay kinematics for both the 1 = 2 and 
1 = channels. This calculation is done on 400 configurations separated by 10 trajectories each. 
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Table 1: Masses of pion and kaons and energies of the two-pion states. Here E' [0 represents the isospin 0, 
two-pion energy when the disconnected graph V is ignored. 



m n E I0 


F' 


Ej2 


m k (0) 


m k (l) 


m k (2) 


0.24267(68) 0.450(17) 


0.4392(59) 


0.5054(15) 


0.4255(6) 


0.5070(6) 


0.6453(7) 



3. Two-pion Scattering 

The 7T — 7r scattering calculation includes 4 contractions, with the name Direct, Cross, Rect- 
angle, and Vacuum diagram [3]. The calculated results from each of these four contractions are 
shown in the left panel of Fig. 1 . Notice that the disconnected (vacuum) graph has an almost con- 
stant eiTor with increasing time separation, so it appears to have an increasing error bar on the log 
plot, while the signal exponentially decreases. The two-pion correlators are fit with a functional 
form Corr(?) = \Z\ 2 {exp{— Et) +exp(—E(T — t)) +C), where the constant comes from the case in 
which the two pions propagate in different time directions. The fitted energies are summarized in 
Tab. 1 . In order to see clearly the effect of the disconnected graph, we also do the calculation for 
the 1 = channel without the disconnected graph. These results are shown with labels that have an 
additional / symbol. The right panel of Fig. 1 shows the resulting effective mass for each case. It 
clearly shows the two pions are attractive in the 1 = channel and repulsive in the 1 = 2 channel. 
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Figure 1: Left: the calculated results for the components of the correlation function defined as Direct(D), 
Cross(C), Rectangle(R), and Vacuum(V). Right: the effective mass for isospin 2 (I2), isospin (To), isospin 
without the disconnected graph (I' ), and twice the pion effective mass (2m K ). 



4. K° to %% Decay Contractions 

The effective hamiltonian for K° to 7T7C decay including the u, d, and s flavors as dynamic 
variables is 

G 10 

H w = -^V^V us Y d Mii) + zy i (n))}Qi (4.1) 
V2 i={ 

where the definition of the ten operators are the same as in [4] , Zi and yi are the Wilson coefficients, 
and T = — V t *V t d/V u dV* s . To calculate the decay amplitudes A2 and Ao, we need to calculate the 
weak matrix elements (7r7r|<2;|^°) on the lattice. 
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For simplicity, we list all possible contractions of (7i7t\Qi\K ^ in Fig. 2. There are 48 differ- 
ent terms, label by circled numbers from 1 to 48, and grouped in terms of geometrical structure into 
typel,type2,type3, mdtype4. The calculation for the correlation functions of (7i7t(t K )\Qi\(t)K° (t^)} 
is then straightforward, for example, in the 1=0 (use notation |/,/ z }) case, 

<00|£ 2 |^°> = /-J={-©- 2 -© + 3 -@ + 3-@-3-@} (4.2) 
v3 

<00|£ 6 |^°> = j\/3{-® + 2-@-© + 2-@ + @-@-©-2-@-@ + @ + @}(4.3) 
A few notes about the contractions shown in the Fig. 2: 

1. The graphs themselves do not carry the minus sign from the odd number of fermion loops. 

2. The dashed line stands for the contraction of colors. If there is no dashed line, it means that 
the trace of color is the same as the trace of spin. 

3. A line stands for a light quark propagator if it is not explicitly labeled with V. 

4. Using Fietz symmetry, it can be shown that there are 12 identities among these contractions, 
such as ©=-©, ©=-©. 

5. Based on charge conjugation symmetry, the average of each of these contractions is real. 

6. The loop contractions in type3 and typeA are calculated with Gaussian stochastic wall sources. 
Two examples of the definition of these graphs: 

© = 7Y{y M (l - Y5)L(x op ,X2)rsL(x2,x p)} • 7Y{y M (l - j 5 )L(x op ,x 1 )y 5 L(x u xo)y 5 S{x ,x op )} 

© = Tr c {Tr s {Yn(l - Y5)L(x op ,x 2 )Y5L(x 2 ,x op )} • Tr s {f{\ - y 5 )L(x p,xi)y 5 L(xi,x )75S(*o,*op)}} 

where xq is the position of the kaon, x\ and x 2 are the position of the two pions, L(x s i n k,x src ) is the 
light quark propagator, and S(x s i n k,x src ) is the strange quark propagator. Tr c stands for color trace, 
Tr s for spin trace, and Tr for both spin and color trace. 

Notice that the type3 and typeA graphs include quark loop integration which results in quadratic 
divergence. However, we should also notice that the operator renormalization allows the mixing 
with the lower dimensional operators s~Ysd and sd [4], where the latter one is forbidden by parity 
conservation. The subtraction of (00\s~Y5d\K } removes the quadratic divergence, and the sub- 
tracted results are calculated as 

(00\Qi\K°) sub = (00\Qi\K°) - on (00\sY 5 d\K°) (4.4) 

(q\q\k°) 

where the subtraction coefficient a,- can be calculated from the K° to vacuum ratio A,. ',,„/, . 

{0\sYsd\K") 

The term (00\sY5d\K ^ comes from two contractions, one connected and one disconnected, 
which are labeled as mix3 and mixA with the coefficent a,- incorporated. To better visualize the 
contributions from different types of contractions, we can write Eq. 4.4 symbolically as 

(QQ\Qi\K Q ) sub = type\+type2 + type3 + typeA — mix3 — mixA 

= type\+ty pel + sub3 + sub A (4.5) 

where sub3 = type3 — mix3, subA = typeA — mixA. 
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Figure 2: All ( x 7i7i\Qi\K ') contractions. They are labelled from left to right, top to bottom by the numbers 1 
to 48. If there is a label '-/+', then it stands for two rows with the '-' sign comes first. If there is an additional 
label '1/s', then it stands for 4 rows with the order '-1', '-s', '+1', '+s', e.g. the first graph is labelled by © 
and ©. 
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Table 2: The fitted results for the weak matrix elements of AI = 1/2 kaon decay. The row label with a prime 
symbol means the disconnected graph is ignored. These are obtained using a source-sink separation of 14. 



i 




1 


2 


3 


4 


5 


Qlatl 


(xl(T 2 ) 


-0.65(38) 


1.75(14) 


1.0(10) 


3.39(80) 


-5.04(91) 


Q\ at ( 


xl(T 2 ) 


-0.4(12) 


1.37(52) 


1.2(33) 


2.9(27) 


-1.7(30) 


i 




6 


7 


8 


9 


10 


Qlatl 


(xl(T 2 ) 


-15.9(10) 


14.35(44) 


44.2(11) 


-1.50(29) 


0.92(29) 


Q lat { 


xl(T 2 ) 


-6.4(40) 


11.6(12) 


34.9(24) 


-1.0(10) 


0.66(97) 



5. Results and Conclusions 

Figure 3 shows the calculated result for operator Q2. The disconnected graph again makes 
huge contribution to the error. We fit the K° to nn correlators with a single free parameter Q\ at {a): 



< n7t(t n )\Qi(t)\K(p) > sub = Q l ?\d)Z* %JZ Z k e- E ^e-^- E ^ 



(5.1) 



where and are calculated from the correlator (K(t)K(O)}, and Z %n and E n% are calculated 
from the two-pion correlator. To see the effect of the disconnected graph, it is ignored and the 
calculated results are shown with an additional / label. The fitted results are shown in Tab. 2. 
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Figure 3: Left: Results of each type contraction of (00(t n = 14)|2 2 (0I^°(^ = °)) 4m (>- Righ t: The tota l 
results of this correlation function and the fitting, where Q2 labels the total result and Q' 7 labels the result 
without the disconnected graph(si<fo4). 

The value of Qi{pi) defined in the MS scheme can be calculated from <2;(jU) = Z^Q^'ia), 
where the Non-Perturbative Renormalization(NPR) matrix Zjj and the choice of the basis follow 
closely those in [4, 5]. The finite volume effect are removed by the Lellouch-Luscher factor (F) 
[6] which relates the quantity M calculated from the lattice in finite volume to the infinite volume 
result A based on our convention of phase space factor: 



|A| 2 = 4n( E ^){p d ^ +<7^}|M| 2 = F 2 |M| 2 



(5.2) 



p J " dp dq 

where p is defined in E nn = \frr^ + p 2 , and q = Lp/2n. Taking the free field limit, this becomes 
|A| 2 = 2(2m K ) 2 mKL 3 \M\ 2 , and the pre-f actors show the different normalization of states in a box 
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and states in infinite volume. For the isospin state in our case, p 2 has a large negative value. 
Even though Eq. 5.2 can be analytically continued to imaginary momentum p, it does not work 
well for large negative p 1 as the function (q) becomes ill defined. We believe that this difficulty 
results because the condition on the interaction range R < L/2 used to derive the Lellouch-Liischer 
factor is not well satisfied. Therefore, it is safer for us to use the free field factor in the 1 = case 
and concentrate on the statistical error from the calculation of the correlators. This problem will 
naturally go away once we explicitly give momentum to the two pions or work in a larger volume. 
Combining everything together, we are ready to calculate the K to 7171 decay amplitudes, 

G 10 

A ie iSl =F^V lld V us ^[(zi(^) + Ty i ^))Z ij (ix)Q\ at (a- 1 )] (5.3) 
V2 i=1 

The calculated Re(Ao) and Im(Ao) are shown in Tab. 3. For comparison, the calculated AI = 3/2 
on shell decay amplitudes of m k = 0.5070(6) to £"72 = 0.5054(15) are Re(A 2 )=5. 394(45) x 10~ 8 
GeV andIm(A 2 )= -0.7792(78) x 10~ 12 GeV. 



Table 3: K° to %% AI = 1 /2 Decay amplitudes in unit of GeV. 



m K 


Re(A> )(xl0-*) 


fe(A )(xl0- 8 ) 


/m(A' )(xlO- 12 ) 


/m(A ) (xl0~ 12 ) 


0.4255(6) 


37.8(2.1) 


28.3(7.8) 


-62.1(5.2) 


-21(20) 


0.5070(6) 


43.5(2.4) 


35.4(9.9) 


-67.7(5.5) 


-48(27) 


on shell 


38.7(2.1) 


30.4(8.5) 


-63.1(5.3) 


-29(22) 



In conclusion, our zero momentum K°C71S MeV) to 71(420 MeV)7r(420 MeV) decay calcu- 
lation gives Re(Ao) with an error of 25%. We find a ratio of Re(Ao)/Re(A2) of roughly 6. Since 
our pion mass is so much heavier than the physical pion mass, a much smaller factor than the ex- 
perimental value of 25 is expected. To have a definite conclusion on Im(Ao) from which we could 
calculate e', we estimate that 4 times statistics could establish a non-zero value if we believe that 
the true result is not too far away from the result without the disconnected graphs. 
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